function plotTP_uvec(u,DU,moonName)

figure;
hold on
[ra,rp]=meshgrid(0:0.0025:1.7);
Tisserand = 2./(ra+rp) + 2*sqrt(2*ra.*rp./(ra+rp));
N = length(u);
for ii = 1:N
    [L1,L2,L3,L4,L5] = librationPoints(u(ii));
    L_ = [L1,L2,L3,L4,L5];
    color = [rand,rand,rand];
    for jj = [1,2,3,4,5]
        % 1. Determine position of lagrange points
        
        L_pos = L_(:,jj);
        
        % 2. determine their jacobi's constants
        
        J_L(jj) = jacobiConst(L_pos,zeros(3,1),u(ii));
        
        % 3. plot using tisserand parameter
        
                contour(ra*DU(ii),rp*DU(ii),real(Tisserand),[J_L,J_L],'Color',color,'LineWidth',3);
%         contour(ra,rp,real(Tisserand),[J_L,J_L],'Color',color,'LineWidth',3);
        
    end
    
    for jj = linspace(J_L(2),J_L(3),7)
        contour(ra,rp,real(Tisserand),[jj,jj],'Color',color,'LineWidth',0.5,'linestyle','-.');
%         contour(ra*DU(ii),rp*DU(ii),real(Tisserand),[jj,jj],'Color',color,'LineWidth',0.5,'linestyle','-.');
    end
%     contour(ra,rp,real(Tisserand),[3.00307873249540,3.00307873249540],'Color','m',...
%                 'LineWidth',3,'linestyle','-.');
% [3.00307873249540]
end
y = linspace(0,1.7*DU(end),10);
% y = linspace(0,1.7,10);

x = @(y) y;
plot(x(y),y)
patch([x(y) fliplr(x(y))], [y max(ylim)*ones(size(y))], 'w')
grid on
for ii = 1:N
    text(DU(ii),DU(ii)+0.3*DU(ii),moonName{ii})
end
xlabel('ra');
ylabel('rp');
set(gca,'FontSize',20)
end